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1 Introduction 

Non self-adjoint eigenvalue problems arise, as is well known, in hydrodynamic and magnetohydrodynamic stability 
[9J [HI [Til HZ] while self-adjoint problems arise mostly in quantum mechanics jl3^. The lack of oscillation theorems in 
the non self-adjoint case makes any computation of the spectrum a very difficult task [lOj . In fact, the eigenvalues 
are scattered over the complex plane and we need first to determine the regions which contain them. A method that 
finds the eigenvalues in a rectangle and in a left half plane has been introduced in [15]. It is based on the argument 
principle with compound matrix method using Magnus expansion. In [1] the authors report on a method that provides 
bounds for the eigenvalues of singular Sturm-Liouville problems over [0, oo) with a complex potential. The method 
consists in obtaining first a floating point approximation to the desired eigenvalue by truncating the infinite interval 
then use interval arithmetic to localize the eigenvalue. In [Q, the author uses the sampling method introduced in [5] 
to compute the eigenvalues of non self-adjoint Sturm-Liouville problems. 

For the mathematical foundation one may consult [HI [TS1 [TS] . On the numerical side [TH1 HH] summarize most of 
the available software dealing with the computation of the eigenvalues of Sturm-Liouville problems. 

In [5], this author introduced the regularized sampling method] a method which is based on Shannon's sampling 
theory but applied to regularized functions. Hence avoiding any (multiple) integration and keeping the number of 
terms in the Cardinal series manageable. It has been demonstrated that the method is capable of delivering higher 
order estimates of the eigenvalues at a very low cost. The purpose in this paper is to extend the domain of application 
of this method to the problem at hand. 



2 Main results 

Consider the following non self-adjoint Sturm-Liouville problem with non-separated parameter dependent boundary 
conditions, 

( -y" + q{x)y = /j, 2 y , i6[0,l] (0 u 

\ A(y(0) 7 y'(0),y(l),y>(l)) T = ^ 

i ,1 , • A/ 1 ( OL\\ (lA 112(m) Cl3(A*) O'liil 1 ) \ i ir, i • i lie,- 

where the matrix A(a) = ; ; ; ; ; ; ; ( has rank 2, and q is a complex-valued function satis- 

V a 2l(W 022(W a23W «24(W / 

fying q € L\ oc (0, 1). We shall not make any assumption on the analyticity of A nor on the growth of its components. 
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The purpose in this paper is to compute the eigenvalues of (|2 . 1 [) with the minimum of effort and a greater precision 
using the newly introduced regularized sampling method [5 , an improvement on the method based on sampling theory 
introduced in [2] . We note here that the analyticity of A and the conditions on the growth of its components imposed 
in [5] are not necessary for the computation of the eigenvalues as shall be seen in the sequel. In fact all what is needed 
is the recovery of certain entire functions hki associated with some base problems defined below. 

It is well known that the spectrum is discrete and scattered over the complex plane which makes difficult its 
computation. Also, there is no result about the distribution nor the multiplicity of the eigenvalues. 

Let y c (x,/j,) and y s (cc,/i), be the solutions of the base problems 

-V" + <l{x)y = i± 2 y , x G [0, 1] . . 

2,(0) = 1 , 2/(0) - [Z - Z) 

and 

-y" + q(x)y = fi 2 y , X G [0, 1] 



y(0) = , 2/(0) = 1 (2 - 3) 
respectively. Then the general solution of the differential equation in (|2.1|) and its derivative are 

y(x, m) = ciy c (x,n) + C2y a (x,fi) 
y'(x,fi) = ciy' c (x,^) + c 2 y' s (x,fj,) 

The boundary condition gives after separating c\ and c 2 , 

c x Awi + c 2 Aw 2 = (2.4) 

where 

wi = (l,0,2/ c (l,Ai),2/c( 1 >A 1 )) T 
w 2 = (0,l,y a (l,n),y'„(l,fj,)) T 

Thus, a necessary and sufficient condition for A = /i 2 to be an eigenvalue is that [i satisfies the characteristic equation 
B(fx) = 0, where B is the characteristic function B{n) = det (Aw\\Aw 2 ) = det[A (wi\w 2 )], that is, 

B(fi) = {a n (fi) +ai 3 (fj,)y c (l,ii) + a u (n)y' c (l, fi)) (a 22 (n) + a 23 (//)2/ s (l,/x) + a 2A (fi)y' s (la, ^)) 
- (a 2 i(n) + a 23 (n)y c (l, /x) + a 24: (fi)y' c (l, fj,)) (a 12 (/x) + a 13 {n)y s (l, /x) + a 14 (n)y' s (l, m)) 

We shall need the following well known results, 

Lemma 2.1 Jlj] sin z/z and cosz are entire as functions of z and satisfy the estimates 

\sinz/z\ </3 e |/m z| /(l + |2|) and |cosz| < e |/m z| 

where (3q = 1.72. 

Using the above lemma one can show the following result to hold. 

Theorem 2.2 y c (x,fj,), y s (x,fj,), y' c (x,/i) and y' a {x,[i) are entire as functions of /i for each fixed x G (0, 1] and 
satisfy the growth conditions, 

\y c {x,ii)-cos(jJ.x)\, \y s (x,fi) -^^\,\y' c (x,fx) + fxsin(fix))\ , \y' s (x, n) - cos(/xx))| < p ie x \ Im »\ 

A* 

for some positive constant Pi . 

In [7] and [5] we have obtained much higher estimates of the eigenvalues than those presented in Theorem 2.2 
above, at the expense of subtracting terms involving multiple integrals. Here and as in [5], we shall stick with the 
estimates given in Theorem 2.2, avoiding any (multiple) integration. We shall show by the same token that we can 
get a higher order estimate of the eigenvalues of the problem at hand at a very low cost. In fact we do not have even 
to keep on increasing the number of sampling points. 
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Let PW a denote the Paley- Wiener space [20] 

PW a = {/entire, |/( M )| < f3e^ lm ^, f \f(^)\ 2 dfi < co} 



Let hki be defined by 



hi2(fi) = 

h.2i(p) = 

Mm) = (^) T "(»i(l,/i)-coBM) 

Then we rewrite y c (l, /u), y' c (l^), y s (l,/i) and y,(l, //) as 



(y c (l,M) - cosm) 
) W (y s (l,M)-^) 
(y' c (l, fi) +fi sin m) 



sin#/£ 
sin 



2/c(1,m) 
2/c(l,M) 



hn(n) 

/l2l(M) 



sin fl/A 
sin 
sin 0/^ 



f COS^i 
, sin 

sin 



^ = MM)(^) m +cos M 



Theorem 2.3 Let $ be a positive constant and m be a positive integer (m > 2). The functions h^i, (fc,Z = l,2) belong 
to the Paley space PW a with a = 1 + and satisfy the estimates 



\hki(fj)\ < 



fh 



\Imu 



k,l = 1,2 /or some positive constant fii. 

Proof: That /ifci are entire and satisfy the given estimates is a direct consequence of Theorem 2.2 and the fact 
that s "g^ M is an entire function of fi and satisfy the estimate in Lemma 2.1^> 

Since the hki(fi) belong to the Paley- Wiener space PW a for each k, I — 1,2, they can be recovered from their 
values at the points fj,j = j—, j G Z, using the following celebrated theorem, 

Theorem 2.4 (Whitaker-Shannon-Kotel'nikov) [20] Let h G PW a , then 

^sinajfi - /ij) 



j = -oo 



Hj = j^. The series converges absolutely and uniformly on compact subsets of C and in L 2 d ^{R). 

For all practical purposes, we consider finite summations, therefore we need to approximate hki by a truncated 
series ftlj . The following lemma gives an estimate for the truncation error. 

Lemma 2.5 (Truncation error) Let h)£p([i) = J2f=-N (Pj) "'g^f-^/ " denote the truncation of hki (m)- Then, 
for \n\ < JVtt/ct, 



hki(p) - h l u ] (p) 



where fo = | j A*"* - 1 ^fcf <A*) I |a - 



< 



J sin/^|/3 3 



^(Tr/o-)" 1 - 1 ^! - 4~ m + 1 



+ 



1 



(iv + \y 



Proof: Since fi m 1 hki{n) 6 L 2 (— oo,oo), Jagerman's result (see [20], Theorem 3.21, p. 90) is applicable and yields 
the given estimate for the hki, k,l = 1, 2^ 

An approximation Bn to the characteristic function B is provided by replacing the hki by its approximation hffl, 
and wc obtain at once, 
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Lemma 2.6 The approximate characteristic function Bjy satisfies the estimate, 



\B(fi) - B N (n)\ < 



sin 9fi 



BpL 



j sin/x|/?4 



7r(7r/<r) m -Vl - 4" m+1 



(iV + 1) 



m-l : 



/or some positive constant ^4. 
We claim the following, 

Theorem 2.7 Let~p 2 be an exact eigenvalue of B of multiplicity n and denote by pL 2 N the corresponding approximation 
of a square of a zero of Bn- Then, for |/ijv| < Ntt/v, we have, 



\fJ-N - Ml 



< 



inf \B( m )(Jj)\ 
1 



sin 9pL 



N 



8fi N 



sin n N |/? 4 



l/n 



7r(7r/<7) m - Vl - 4- m+1 
l/n 



1 



1 



where the inf is taken over a ball centered at /xjv with radius — Jt| and noi containing a multiple ofir/9. 
Proof Since /I is a zero of -B with multiplicity n, then 

b(m)-^)= (P ~^ ) " ^")(m) 



for some /i. Thus, 



to!|£(7I)- 5(^)1 



|a(»0(/Z)| 



< 



inf |-B( m )(/2)| 
1 



sin #/i 



A 



(9// 



A 



sin^Ar|/3 4 



7r(7r/CT) m - 1 \/l - 4~ m 



+1 



(7V + l)m-l 



where the inf is taken over a ball centered at /ijv with radius |/ijv — /i| and not containing a multiple of 7r/#. Thus, 
the result 



3 Numerical examples 

In this section, we shall present a few examples to illustrate our method. We have taken = 1/(N — m) in order 
to avoid the first singularity of ( s "g^ w ^ ■ The sampling values were obtained using the Fehlberg 4-5 order Runge- 

Kutta method. The first two problems are taken from j¥j in which the authors use interval arithmetic to localize the 
eigenvalues of singular Sturm-Liouville problems with complex potentials. The third problem, taken from pQ, shows 
that the regularized sampling method provides much better results than the sampling method without regularization. 
The last example demonstrates that our method can estimate the eigenvalues with a great precision even in situation 
where other methods might introduce spurious eigenvalues and/or miss some of them. We shall mention however that 
we shall not make use of the error estimate given above for the time being. The method consists first in the recovery 
of the entire functions h^i with great precision, then use the boundary conditions to determine the characteristic 
function. The zeros of this characteristic function are the square roots of the sought eigenvalues. We shall denote 

Example 3.1 (Taken from [J]) Consider the singular Sturm-Liouville problem 

{—y"(x) + 10 Lsn\xe~ x y(x) = Xy{x) , < x < 00, 
1/(0) = 
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We shall use interval truncation and compute the eigenvalues of 

J —y"(x) + 10 Lsuvxe~ x y~ l (x) = /i 2 y 7 (x) , < x < 7, 

I 2/7 (°) = ' ^(7) = WAyW 

and as in [4] we shall take 7 = 10. The second boundary condition has been obtained by considering the Jost solution 
y = e L ^ x and its derivative y' = i/ie tAia: , thus, 2/(7) = t/x 3/(7). In [4] the authors obtained an eigenvalue lying in 
1.6043912ft + 1.7978849^ 1 where the notation 2.l| stands for the interval [2.14,2.16]. Taking N = 40, and for 
different values of m we obtained the results summarized in Table 1, 



m 


Approximate Eigenvalue 


5 


1.604391348283 + 1.797884747658 1 


10 


1.604391251270 + 1.797884973775 1 


15 


1.604391251323 + 1.797884973746 1 



Table 1: Approximation of an eigenvalue for different values of m in Example 3.1 

Example 3. 2 (Taken from [4]) Consider the singular problem 

J —y"(x) + 10 l e~ x y(x) = \y(x) , < x < 00, 
\ 2/(0) = 

We shall use interval truncation and compute the eigenvalues of 

f -y"(x) + 10 ie~ x Vl {x) = [P-y^ix] , < x < 7, 
1 2/7(0) = . 2/7(7) = ^2/7(7) 

and as in [4] we shall take 7 = 10. In [4] the authors obtained an eigenvalue lying in 2.8122672g| + 2. 17223818^ l. 
Taking N — 40, and for different values of m we obtained the results summarized in Table 2, 



m 


Approximate Eigenvalue 


5 


2.812264032443898167911 + 2.17223723666731852353 1 


10 


2.8122672894628469454261 + 2.172238191264223861 1 


15 


2.812267288417814133626 + 2.172238191179093864 t 



Table 2: Approximation of an eigenvalue for different values of m in Example 3.2 

Example 3.3(Taken from p]) Consider the non self-adjoint problem 

f -y"(x) + (3 - 2 L)y(x) = Xy(x) ,0<x<ir, 
\ 2/(0) = y(n) = 

The exact eigenvalues of the original problem are Xk = k 2 + 3 — 2 t, k = 1,2, .... Taking N — 40 and m = 10, we 
obtained the results summarized in Table 3, 
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15 
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1.05xl0~ 9 


16 


259-2 i 


258.999999996079-1.9999999963552374114 i 


5.35xl0~ y 


17 


292-2 t 


291.999999991669-1.9999999961982129591 l 


9.15xlO~ M 


18 


327-2 l 


327.000000004474-2.0000000150098733387 i 


1.56x10-* 


19 


364-2 t 


364.000000082618-2.0000000798746168342 i 


1.14x10"' 


20 


403-2 i 


403.000000232131-2.0000001291033812382 i 


2.65xl0~ 7 



Table 3: Exact and Approximate eigenvalues in Example 3.3 



Example 3.4 Consider now the following non self-adjoint Sturm-Liouville problem with complex potential and 
parameter dependent boundary condition, 

-y"{x) + e 2 LX y(x) = Li 2 y{x) , < x < 1, 
y(0) + M y(lj = 
y'(0) = 

Here again we are in a position to derive the exact characteristic function which in fact can be expressed in terms of 
Bessel functions. Indeed, let A = /j 2 and consider the change of variables t — e LX . The differential equation becomes 
the Bessel equation of order /i given by 

whose solution is 

z(t) = aJpit) + aJ-nit), 

where and J7_ M are the Bessel functions of the first kind of order fi. 
Returning to the original variables, we obtain 

y(x) = c 1 J fl (e' x ) + c 2 J-^ x ). 

Taking into account the boundary conditions, we obtain the homogeneous system in ciand c 2 

f ci + c 2 J-^(l) + fi (cM&) + c 2 J-^&)) = 

1 Cl j;(i) + C2 j^(i) = o 

In order to have a nontrivial solution, a necessary and sufficient condition is to have B exact {[i) — where 
is the characteristic function. Now, using the well known result 

4:M X ) = (^-m-iW - /2, 
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we obtain 

Bexact(v) = det 

Taking N = 40, and m = 10 we obtained the results summarized in Table 4. 



J M (l)+/i^(e t ) J_, I (l) + M^- Al (e t ) 

(^-i(l)-^+i(l))/2 (J M (l)-J- M+ i(l))/2 



Approximate Eigenvalue 



Absolute Error 



Index 



Exact Eigenvalue 



5.549xl0" 1:> 
3.393X10" 14 
3.977xl0" 13 
8.004x10" 
2.064x10 



Relative Error 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



4.9685430929323576+0.3906545895360696 t 
20.60271034889337+0.75023252353154 i 
64.14038244804547+0.68422837531133 i 
119.34792168887388+0.71497240479401 t 
202.31443747778734+0.70057212586525 t 
419.44558800598641+0.70446189520144 i 
553.61789373762934+0.70954623577257 i 
715.53365857906959+0.70595783818772 i 
889.18520034251622+0.70898948206981 i 
1090.57859485902126+0.70668585309098 I 
1303.70898166607058+0.70869788000992 i 
1544.58037965386611+0.70709389168016 I 
1797.18943505543540+0.70852627026801 f. 
2077.53900632820814+0.70734525957323 o 
2369.6266391592291618+0.70841680475450 
2689.45447190894724+0.70751097714777 I 
3021.02063035583927+0.70834272705249 i 
3380.32677490847313+0.70762595333135 I 
3751.37142735725052+0.70829027475581 I 
4150.15591451714336+0.70770896968510 t 
4560.67904058883973+0.70825178027457 o 
4998.94189026423779+0.70777085938814 i 
5448.94347623327640+0.70822269647649 i 
5926.68470186115217+0.70781822661900 i 
6416.16473814590617+0.70820018792052 i 



4.9685430929323625+0.3906545895360721 
20.60271034889340+0.75023252353155 i 
64.14038244804526+0.68422837531099 t 
119.34792168887345+0.71497240479334 i 
202.31443747778739+0.70057212586545 I 
419.44558800598892+0.70446189520528 I 
553.61789373762976+0.70954623577282 i 
715.53365857906140+0.70595783817453 t 
889.18520034250143+0.70898948204681 I 
1090.57859485902214+0.70668585309385 t 
1303.70898166611992+0.70869788008925 i 
1544.58037965396658+0.70709389183628 i 
1797. 18943505544546+0.70852627027458 i 
2077.53900632774381+0.70734525883073 i 
2369.62663915816209+0.70841680308871 i 
2689.45447190899851+0.70751097732396 t 
3021.02063036245257+0.70834273761505 i 
3380.32677492777072+0.70762598345876 ( 
3751.37142735927201+0.70829027606949 ( 
4150.15591430123958+0.70770862603678 ( 
4560.67903973138474+0.70825043937887 ( 
4998.94189032279592+0.70777105963898 i 
5448.94349859286012+0.70825825178573 ( 
5926.68487236793343+0.70808488449668 ( 
6416.16478405538947+0.70823000684001 ( 



1-13 
-13 



4.582X10" 12 
4.969X10" 13 
1.553X10" 11 
2.734X10" 11 
3.00X10" 12 
9.341X10" 11 
1.856X10" 10 
1.201X10" 11 
8.757xl0" 10 
1.978X10" 9 
1.834xl0" 10 
1.246xl0" 8 
3.577xl0" 8 
2.410x10"° 
4.058X10" 7 
1.591X10" 6 
2.086X10" 7 
0.00004200 
0.0003165 
0.00005474 



1.113x10" 
1.645x10" 
6.201x10" 
6.706x10" 
1.020x10" 
1.092x10" 
8.977x10" 
2.170x10" 
3.075x10" 
2.751x10" 
7.165x10" 
1.202x10" 
6.687x10" 
4.215x10" 
8.348x10" 
6.822x10" 
4.125x10" 
1.058x10" 
6.426x10" 
9.778x10" 
3.489x10" 
4.173x10" 
7.708x10" 
5.340x10" 
8.532x10" 



Table 4: Exact and Approximate eigenvalues in Example 3.4 



4 Conclusion 

In this paper, we have used the regularized sampling method introduced recently [5] to compute the eigenvalues of non 
self-adjoint Sturm-Liouvillc problems with nonseparable parameter dependent boundary conditions. We recall that 
this method constitutes an improvement upon the method based on Shannon's sampling theory introduced in [2] since 
it uses a regularization avoiding any multiple integration. The method allows us to get higher order estimates of the 
eigenvalues at a very low cost. We have presented a few examples, including singular ones, to illustrate the method 
and compared the computed eigenvalues with the exact ones when they are available. 
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